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Abstract 

We consider a highly anisotropic d = 2 Ising spin model whose precise definition 
can be found at the beginning of Section[2] In this model the spins on a same horizontal 
line (layer) interact via a d = 1 Kac potential while the vertical interaction is between 
nearest neighbors, both interactions being ferromagnetic. The temperature is set equal 
to 1 which is the mean field critical value, so that the mean field liirrit for the Kac 
potential alone does not have a spontaneous magnetization. We compute the phase 
diagram of the full system in the Lebowitz-Penrose limit showing that due to the 
vertical interaction it has a spontaneous magnetization. The result is not covered 
by the Lebowitz-Penrose theory because our Kac potential has support on regions of 
positive codimension. 


1 Introduction 

This work focuses on the proof of the Lebowitz-Penrose limit for a highly anisotropic 
d = 2 Ising spin model which has been first studied in [2], its precise definition can be 
found at the beginning of Section [2j In this model the spins on a same horizontal line 
(layer) interact via a d = 1 Kac potential while the vertical interaction is between nearest 
neighbors, both interactions are ferromagnetic. The temperature is set equal to 1 which 
is the mean field critical value (without vertical interactions), so that the mean field limit 
for the Kac potential alone does not have spontaneous magnetization. However in [2] it is 
proved that even a small vertical interaction is sufficient to produce a phase transition at 
least for small values of the Kac scaling parameter 7 . The idea in [2] is to study a model 
with fewer vertical interactions (those left have a chessboard structure): by the Ginibre 
inequalities if a spontaneous magnetization is present in the reduced model then it is also 
present in the true system as well. The advantage of working in the reduced system is 
that one can split the system into blocks of two layers, the vertical interaction is left only 
inside each block so that blocks do not interact vertically with each other; the horizontal 
interaction is unchanged. As a consequence in [ 2 ] it is shown that it is sufficient to carry 
out the Lebowitz-Penrose coarse graining procedure only for two-layer systems. It is then 
proved that this can be done, that there is a positive spontaneous magnetization in the 
limits volume to infinity and then 7 —)• 0 and that such a property remains valid also at 
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finite small 7 > 0 . However the value of the spontaneous magnetization for the reduced 
system is certainly smaller than the real one because half of the vertical interactions has 
been dropped. 

The problem of studying directly the original system and in particular to hnd its true 
spontaneous magnetization has been left open in [2], we attack it in this paper determining 
explicitly the limit phase diagram of the true system when first the volume goes to infinity 
and then 7 —)• 0 . This is not covered by the Lebowitz-Penrose theory because our Kac 
potential is singular having support on regions of positive codimension. We hope in a 
successive paper to prove that there is a positive spontaneous magnetization also at 7 > 0 
which converges as 7 0 to the one found here. 

This work is part of a more general project (which besides us involves several other col¬ 
leagues) where we want to study systems with Kac potentials having support on regions 
of positive codimension plus short range interactions, both in equilibrium and non equilib¬ 
rium. The description of the system is hybrid: referring to our Ising model we can make 
a coarse graining on each layer and introduce macroscopic variables but the interaction 
between layers is microscopic and it is described by an effective interaction to be deter¬ 
mined. The purpose is to derive such an effective hamiltonian and find its ground states. 
In this paper we compute the limit ground state energy but we hope in the future to study 
the excited states and derive the large deviations functional. 

Similar structures are present in SOS models, for instance in the SOS interface models 
where the real valued spin variables Sx, x G represent the position of an interface 
in Evidently the model is obtained by an anisotropic scaling limit for which the 

interface becomes sharp,the point Sx on the “vertical” line through x, while the interaction 
among spins remains short range. We hope to establish such connections starting from 
models like the one considered here. 

The appearance of a macroscopic description on the layers may also originate from a 
canonical constraint with or without the presence of a Kac potential. Considering the 
system in a hnite box we may fix on each horizontal layer the total magnetization which 
gives rise to a multi-canonical ensemble. Indeed when we study the system with Kac 
potentials following the Lebowitz-Penrose procedure we coarse grain and get such multi- 
canonical ensembles. Our analysis will be based on a proof that equivalence of ensembles 
extends to such cases. 

The multi-canonical constraint appears naturally in dynamical problems when we consider 
a Kawasaki dynamics on each layer so that the total magnetization on each layer is constant 
in time. The vertical interaction affects the rates of horizontal exchanges on the layers 
so that in the hydrodynamic limit the evolution is conjectured to be ruled by coupled 
diffusions. An interesting variant would be a weakly asymmetric simple exclusion on each 
layer with small interactions among layers which should be in the KPZ class of systems. 
We refer to the introductions in m-m for more references and a list of open problems and 
conjectures, in particular the connection with quantum Ising models (via Feynman-Kac), 
phase transitions for the hard-rods Kac-Helfand model and the dependence on 7 of the 
critical value of the vertical interaction for a phase transition to occur. 

We conclude the introduction by observing that highly anisotropic interactions are present 
in nature, the best example is the graphite where horizontal structures are rather free to 
slide one with respect to the other. However it may happen that even a small interaction 
among layers produces macroscopic effects. For instance for bilayer graphene samples 
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interacting via an interlayer coupling constant mum the presence of a band gap in the 
energy spectrum, which is tunable by an external electric field, paves the way to a variety 
of applications in electronics [8]. 

Multilayer graphene samples have also gone, recently, under intense investigation [3 do], 
which revealed the rise of exceptional thermal conduction properties for these materials 
as well as the possibility to control the thermodynamically stable cristalline structure of 
the material through an external voltage. 


2 The model and the main result 


As mentioned in the introduction one of our aims is the extension of the Lebowitz-Penrose 
theory to cases where the support of the Kac potential has a positive codimension. This 
is what we accomplish in this paper in the simple context of the d = 2 Ising model. Let A 
be a square in 7?, L its side, {x,i) its points. Write a G {—1,1}^ for a spin configuration 
in A, define a{x, L + 1) = a{x, 1), a{x + L,i) = a{x, i) and let 


KL.Li-) 

= + Hr\a) + 




2=1 


L L 

Hr\a) 

= y{-Ay (T(x,i)(T(x,z + 1)} 

ai=l 2=1 

^/lextA(f^) 

= - y /iext<7(x,i) 


(x,i)eA 


^^ 7 , 1 , (o') is the Kac hamiltonian, it has only horizontal interactions; is the hamil- 

tonian of a nearest neighbor Ising model with only vertical interactions; L^hext,i(o') is the 
energy due to the external magnetic field /iext ■ We suppose that 


J'y{x,y) = c^-fj{-fx,-fy) 


( 2 . 2 ) 


where J(r, r') is a smooth, symmetric probability kernel on M which vanishes for |r—r'| > 1; 
is such that 


J2Jj{x,y) = 1 


(2.3) 


Since f J{r,r')dr' 


1, c.y —)• 1 as 7 —)■ 0. Let 


■y,hext,L / ^ 


(o') 


(2.4) 


be the partition function relative to the hamiltonian ^ i,(o')- Call f\{m) the free 

energy density with magnetization density m relative to the hamiltonian since 

the horizontal interactions are absent f\{m) is equal to the free energy of the d = 1 Ising 
model with only nearest neighbor interactions of strength A. 
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Theorem 1 For A > 0 small enough 

log ZP'f r r m 2 -1 

lim lim -= - inf | -/lextm + [—^ + /aM]| (2.5) 

7^0 L—>-oo |A| mG[—1,1] 2 J 

After a few comments on Theorem [1] we give a heuristic derivation of (j2.5p followed by a 
description of how proofs are organized in the various sections. 


2.1 Remarks on Theorem [T]. 


• (12.5p remains valid for general Van Hove regions and boundary conditions since the 
interaction has hnite range for any fixed value of 7 > 0. The restriction to small A 
is needed for cluster expansion, it is technical and could be presumably removed. 

• The limit in (12.5p is the sum of the external magnetic field energy —/lextm, the mean 
field energy —m?/2 and the vertical free energy fx{m): it reflects the analogous 
splitting of the hamiltonian in ( 12 .ip . 

• limi^oo g 77 jext,-L _. ig pressure of the system with hamiltonian L 

By ferromagnetic inequalities P'y{hext) is for any 7 > 0 a convex function of hext dif¬ 
ferentiable at any hext 7 ^ 0 ; its derivative is the magnetization which is equal to the 
average spin for the unique DLR measure at the given values of hext and 7 . The lim¬ 
its (by subsequences) of P^{hext) as 7 —)• 0 are thus convex functions and Theorem [T] 
proves that the limit actually exists (without going to subsequences) and identifies 
its value. 


• The limit of P^{hext) as 7 —>• 0 is the pressure P{hext) in the Lebowitz-Penrose 

limit when first |A| —)• 00 and then 7 —)• 0 . ()2.5p shows that P{hext) is the Legendre 

2 

transform of the function [—^Pfx{m)] and therefore the free energy Fx{m) defined 
as the Legendre transform of the pressure P(hext) is equal to the convex envelope: 

2 

Fx{m) = Cnj - ^ + fx{m)] (2.6) 


• (12.61) is in agreement with the Lebowitz-Penrose result which states that the limit 

2 

free energy density is the convex envelope of — ^ plus the free energy density of 
the reference system (i.e. without the Kac interaction). The Lebowitz-Penrose anal¬ 
ysis however applies if the Kac interaction is non degenerate being positive in two 
dimensional regions. Our theorem shows that this is not necessary. 

• When A = 0, /o(m) = —S{m) where S{m) is the entropy of the free Ising model: 

, 1 -l-m, 1 -l-m 1 -m, 1 -m 

-5(^) = __log__ + __log^— (2.7) 

2 

In this case —^ + fo{m) is strictly convex and coincides with Foim). When A > 0 

2 

we shall see that the function —h fx{m) is no longer convex. In fact the Taylor 

expansion of —S{m) gives 


— S{m) 


log 2 -F ^ 
k=0 


—A - ^ 

2k + l2k + 2 


( 2 . 8 ) 
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and to leading orders in A, f\{m) = —S{m) — \m? so that —^ fx{fn) has a 

double well shape with minima at ±\/^ and F\{m) is constant in the interval with 
endpoints ±\/6A. The spontaneous magnetization is then \/6A to be compared with 
the value \/^ found in [2] for the system with reduced vertical interactions, as 
described in the introduction. 

• The proof of Theorem [T] does not require the use of a non local free energy functional 
as the one introduced by Lebowitz-Penrose, but we have nonetheless established some 
basic ingredients for its derivation which will be used in a future work to study the 
large deviations. 


2.2 Heuristic derivation of the mean field equation 

Let {a{x,i)) =: m be the average spin in an extremal, translation invariant DLR measure 
at 7 > 0. Then 

((t(x, i)) = (tanh{y~] J^{x, y)a{y, i) + A[cr(x, i + 1) + a{x, i - 1)] + hext}) (2.9) 
y 

By the law of large numbers Yhy J-yix, y){cr{y, z) — m) —)• 0 in the limit 7 —)• 0, recall that 
J-y{x,y) = 1. In such an approximation (12.9p becomes 

{a{x, i)) = (tanh{A[cr(x, i + 1) + a{x, i - 1)] + hext + rn}) (2-10) 

This is the equation for the average spin in a d = 1 Ising model with only nearest neighbor 
interactions of strength A and magnetic field hext + Then the average spin is equal 
to the thermodynamic magnetization m which is related to the free energy f\{m) by a 
variational principle which gives 

d / TTX^ \ 

0 = (hext + m)- fx{m) = hext “ ^ j ( 2 . 11 ) 

in agreement with (j2.5p - (j2.6p . 

2.3 Organization of the paper 

The proof of (12.5|) distinguishes large and small values of the magnetization and conse¬ 
quently of the magnetic field. Large magnetic fields are studied in Section [3] by using the 
Dobrushin high temperature theory based on the Vaserstein distance; the “small” values 
of the magnetic field are studied in the remaining sections. In Section U] we give the scheme 
of proof of Theorem [T] which is based on the following steps (each step being discussed in 
a subsection). (1) a coarse graining procedure a la Lebowitz-Penrose to reduce to a d = 1 
system with only nearest neighbor interactions and without Kac potentials. The price is 
that we have a variational problem with multiple constraints as we have fixed the mag¬ 
netization on each layer. (2) We then consider the analogous problem in the multi gran 
canonical ensemble where on each layer we have a magnetic field. The partition function 
of such a system is studied in details using cluster expansion under the assumption that 
A is sufficiently small. (3) We prove an extended equivalence of ensembles so that the 
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original variational problem with constraints given by the magnetization is replaced by a 
variational problem where one needs to optimize on the value of the auxiliary magnetic 
fields. (4) The proof proceeds by showing that the minimizer is made by magnetic fields 
equal to each other on each layer. (5) We then show that Theorem [1] follows. 

In Section [5] we prove a combinatorial lemma which says that any monomial in 

the variables ui, ..,Uk, ni + ■ ■ ■ Uk = N > 2 , can be written as a sum of one body mono¬ 
mials Piuf, Pi positive numbers, plus a sum of terms proportional to gradients squared, 
'^i<j dij polynomials of degree N — 2. This is the essential property 

needed to prove that the minimizers are homogeneous. 

The proofs of all the above statements are reported in successive appendices. 


3 Large magnetic fields 

The heuristic argument presented in Subsection 12.21 is made rigorous for large magnetic 
fields in the following theorem. 


Theorem 2 For any X > 0 let hext > 0 be so large that 

1-F2A 1 

T - <C — 

cosh^(/iext — 1 — 2A) 4 


(3.1) 


Then (i) for any 7 > 0 there is a unique DLR measure (by ferromagnetic inequalities the 
statement actually holds for any hext 7 ^ Oj; (ii) its magnetization (the average value 
of a spin) converges as 7 —)• 0 to the value m for which ( 12 .lip holds; (Hi) m is the unique 
minimizer of (I23|) and 


lim lim 

7^0 L—>-cx) 


log^rs. 


7!^ext iL 


|A| 


2 

'1 

- hextm -g [—— + f\im)]j 


(3.2) 


As we have already mentioned the proof of Theorem [2] is based on the techniques intro¬ 
duced by Dobrushin to prove his famous large temperature uniqueness theorem. In this 
way we get uniqueness of the DLR measures and exponential decay of correlations for any 
fixed 7 > 0. We then use an interpolation procedure to derive the phase diagram of the 
system which was introduced in [3] to study the corrections in 7 to the mean field limit 
and thus prove Theorem [2j The details are reported in Appendix lAl 


4 Theorem [Tt scheme of proof 


Theorem [T] is thus proved for large magnetic fields and the remaining part of the paper 
deals with the “bounded” magnetic fields. To be precise we suppose hereafter A G (0,1), 
but further requests on the smallness of A will be asked later on, and restrict to magnetic 
fields 


hext e [0,h*], 


h* := 


_ 

cosh^(/i* 


3) 


1 

4 


(4.1) 


6 







as Theorem [2] covers the values hg^t > h*. By default in the sequel /lext £ [0, h*] (the 
analysis of negative magnetic follows by symmetry). 

The first step is to use coarse graining as in Lebowitz-Penrose. 


4.1 The Lebowitz-Penrose procedure 

In this subsection we use the Lebowitz-Penrose procedure to reduce to a partition function 
where the Kac potential is absent. Let us first recall the Lebowitz-Penrose result and 
consider the partition function ^ ^ with the same short range, vertical interaction as 
in our case (the “reference system” in the Lebowitz-Penrose terminology) but with a Kac 
potential which has support on a region of full dimension {d = 2 ). After coarse graining 
and exploiting (i) the smoothness of the Kac potential, (ii) the ferromagnetic nature of the 
Kac potential, Lebowitz-Penrose have proved that Z^^^ ^ ^ has the same “asymptotics” 
as 



where A is a square of side i, i the integer part of 7 and 


-^A — {-1, -1 + . . . , 1 - -r^, 1} 


the set of all possible values of the empirical spin magnetization in A. 
By same “asymptotics” we mean that 



(4.2) 


The same procedure works in our case as well leading to Theorem [3] below whose proof is 
given in Appendix iBl 


Theorem 3 Let A and i be as above, = {—1, —1 , 1 — 1}, m^{x, i), (x, i) G 

A, a function with values in whieh depends only on i, 



(4.3) 


CTI 




Then there is m+ G (0,1) so that 


lim lim -—- log Z 

7^0L^oo |A| 



lim log Z\ 


A 


■max 


(4.4) 


where 


\ogZ^ 


max 


+heyitm{x,i) - (4.5) 


a:,2)GA 
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(j4.2l) and (|4.4p are identical but the meaning of is different in the two cases. In (I4.2p 
it is a max over a scalar m of the canonical partition function with magnetization m. By 
classical results on the thermodynamic limit this is related to the free energy of the system 
and one gets a formula as on the right hand side of (12.5p . Thus one has essentially finished 
once he gets ()4.2p . in our case instead (14.4p is just the beginning of the work. In fact the 
variational problem behind (|4.5p involves a vector ruA in a space whose dimensions diverge 
in the thermodynamic limit. Moreover the relation between and the d = 1 free 

energy f\{m) which appears in (12.51) is not evident due to the multi-canonical constraint 
of fixing the magnetization on each layer. 

The picture looks simpler if we replace the multi-canonical ensemble by a gran canonical 
ensemble with auxiliary magnetic fields on each layer: let then h = (hi,.., hi) and 



(4.6) 




The goal is to rewrite (j)i{mj\) in terms logZ^^ h thus proving an extended version of 
the equivalence of ensembles theorem. The first step in this direction is to get a full 
understanding of Za,/i as provided by the cluster expansion. 

4.2 Cluster expansion 

We first observe that 


log ZA,h = i log Zi^h 


^i,h= a{£ + l) = ail) (4.7) 




with Zi^h the partition function of the d = 1 Ising model with nearest neighbor interactions 
of strength A and space dependent magnetic field h. We define 


£ 



(4.8) 


2=1 

Ui := tanh{hj} 


(4.9) 


In Appendix Owe shall suppose A small and use cluster expansion to prove: 


Theorem 4 For any A > 0 small enough 



(4.10) 


7V(.) 


where N{-) : [1,1] —)• N and 


i 



(4.11) 


i=l 


The coefficients Aiv(.) satisfy the following bounds. Call 


:= A-5A2, \N{.)\ =Y,N{x), ||iV(.)|| = max{|iV(.)|, i2(iV(-))} 


X 


(4.12) 


where R{N{-)) denotes the eardinality of the support ofN{-) (i.e. the smallest interval 
whieh contains the set {i : N{i) > 0}J. Then for any i € [1,-^] and any positive integer M 


E (4.13) 

Ar(-):Af(i)>0,||Af(-)||>M 

Moreover ^Ar(.) = 0 if |A^(-)| is odd and there are eoeffieients ak, k > 0, and c so that 

^ (4.14) 

N(-):\N{-)\=2 i<j 

|ai — A| < cAe“^^, \aj-i\ < (4-15) 


4.3 Equivalence of ensembles 

The magnetizations associated to Zi^h, as defined in (I4.7p . are m = (mi, ..,m£) 

mi = Zi^h (4.16) 

Ohi 

which are thus expressed via h in terms of (ui, We write more explicitly (I4.16p as 

mi = Ui + ^i{u), i{u) = {I - uj) ^ (4.17) 

N(-y.N{i)>0 

with N^^\k) = N{k) for /c / f and N^^\i) = N{i) — 1. In Appendix |F] we will prove 
that there is a one to one correspondence between u and m so that we may write u as a 
function of m. 

Theorem 5 For any A > 0 small enough the following holds. For any m such that 
\mi\ ^ " 1 + (fn+ as in (|4.5p ) there is a unique h such that (14.161) holds for any i = 1, 
and there exists h^ > 0 so that all the eomponents of h are bounded by h^. 


Theorem 6 For any A > 0 small enough the following holds. For any = m = 
(mi,..,m^), \mi\ < m+, i = callh = {hi,..,h£) the magnetic fields associated tom 

via Theorem\^ then for any a G (^, 1) there is c so that 


^ log{e Za, ft} + 0 £(mA) 


< c£^-^ 


(4.18) 


where (j)e{mj\) is defined in (|4.3I) and Z^ ,ft in (jM]). 


As a consequence: 
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Theorem 7 Let he as in dMD, then 


lim 
| A |—>-00 


logZ^- 


lim — max 

£—>■00 l }v.\hi\<h^ 


+ (/^ex. 


hi)mi + log 


(4.19) 


where mi is the function of h defined in (I4.16p - (l4.17j) . 


By ()4.4p and (I4.19p we get 


1 


1 


lim lim 7—7 log r = lini — max 

7^0L^oo |A| 7Aext,L £^00 £ h 


+ (^ext - hi)mi + log Ze^h] (4.20) 


i=l 


4.4 The quadratic structure of the effective hamiltonian 

The goal therefore is to study the ground state energy of the effective hamiltonian 
£ 2 

^£% ~ ~ — himi + hext^Tii} — log(e^* + e ^*) — log Zf^^ + Aq (4-21) 

i=l 

regarded as a function of n = (tti, For convenience in (I4.2ip we have subtracted to 

log Zfi^ (which is defined in (jd.lOp ) the first term of the expansion (Id.lOj) (with N{-) = 0), 
which is a constant. 

By Theorem[5]we can restrict to the set of u : \ui\ < = tanh(/i+),i = 1,and in the 

sequel we will tacitly restrict to such a set. We will prove that the inf over u of is 
achieved by vectors u with all components equal to each other which is maybe the most 
relevant/original result of this paper. 

We start by making explicit the leading terms in (|4.2ip for A small. To this end and 
recalling that log(e^* +e“^*) = hiUi + S{ui), the entropy S{u) being defined in (j2.7[) - ()2.8p . 
we write 


log(e^' + e ^') = hiUi + ^ + T{ui) 

T{u)= log2 + y / / ^2^+2 

^ ^ ^ ^ 2k + 12k+ 2 

(4.22) 

k=l 

X ^ 

log -Ao = -- ^{ui+i - Ui)"^ + 0 

i=l 

ii := {hi - Ui){l - uf) 

(4.23) 

(4.24) 

4'i = A(1 - uj){ui+i + Ui-i) + 

(4.25) 


where 4'* is defined in (14.171) and 0 and 4** are defined by the above equations. By some 
simple algebra, see Appendix [G] for details, we can rewrite as: 
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Lemma 1 With the above notation: 
i 

^1% = ^ {T{Ui) - hextUi - 2Xhext[Ui - uf] + 2X^iUi} 
i=l 

^ ( X 'I'? I 

“1“ ^ ^ ^ '^i+l') “1“ {hi hext^i r © 

i=l ^ ^ 

e 

^ ^ {hext{ui -|- {ui-\-l Mj) + ) (?/); (4.26) 

i=l 

The terms with 0, T? and are “under control” in the following sense: 

Theorem 8 Call 

(4.27) 

Then for any h^xt S [0, h*] there is a continuous function f^^\u) on [—1,1] and functions 
b\^j{u), i < j, so that 

i 

= Y1 + Y1 

i=l l<i<j<e 

with 

i 

^ bf]{u){ui - Ujf > -cA^+i ^(ui - Ui+if (4.29) 

Wi<3<t- *=1 

and c a positive constant. 


The proof of Theorem [ 8 ] starts from (|4.10p and it is based on a representation of the 
monomials u^b) as sum of one body and gradients squared terms which is established in 
Section [5l After that we exploit the properties of the coefficients ^Af(-) stated in (14.1,ip . 
(I4.14P and (I4.15p . The computations are straightforward but lengthy, the details are 
reported in Appendix [Hi 

Define 

ei{u) := ~ (4.30) 

^ 2+1 

when Ui ^ ttj+i and equal to df^/du{v) when Ui = Ui+i = v. Then 6 i{u) depends only 
on Ui and Uj+i, and as a function of Ui and Uj+i is continuous, symmetric and bounded 
in |uj| < R, |ui+i| < R, R < 1. It is essential in our proof that 9i{u) < We checked 
numerically that this is “fortunately” true and indeed we found a rigorous proof, reported 
in Appendix m of an upper bound smaller than 1/2: 
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Proposition 1 


< g 


(4.31) 


It follows from (j4.26p and (|4.28p that 

^ i 


treff _ 

— 


^ /(Wi) + ^ X] I ^ “ ^ext(Wi + Ui+l) - 9i{u) \ {Ui - Ui+if 
i=l i=l ^ ^ 


+ ^2 %! - Ujf 

i<j 


where 


(4.32) 


f{ui) = {T{ui) - hextUi - 2XheKt[ui - uf] + 2X^iUi} + (4.33) 

Let u* be the minimizer of /(•); then by (I4.32p 

inf < I inf f{u) =: if{u*) 


as the right hand side is the value obtained by choosing all Ui = u*. 
Let 




1 3 

2 ~ 8 


then if hext £ [0, ho] using (14.291) we get 


2^12 

i=l 2=1 


(4.34) 






2^12 

2=1 2=1 ^ 


(4.35) 


Hence 


infi2?i>^/(u*) 


for A so small that 2cA3 < 2 — because by (14.351) we then get a lower bound by neglecting 
the sum of the terms with (li, — lij+i)^. We have thus proved: 


efr 


Theorem 9 Let hext G [0, Lq] and X be so small that 2 cA3 < 2 ~ then the inf of 
is equal to the min of over homogeneous h, namely h with all its components equal 
to each other. 


Thus the inf of is achieved when all the components of h are equal to each other, in 
such a case we can compute explicitly the minimizer, see the next subsection. The result 
comes from the quadratic structure of the hamiltonian, (I4.28I) - (I4.29D . somehow reminiscent 
of the Ginzburg-Landau functional whose integrand has the form W{u) + |Vrtp and the 
gradient term forces the minimizer to be a constant. 
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The argument used to prove Theorem [S] does not extend to the complementary case when 
hext € [ho, h*] because {^ — hext{ui + rti+i)} in (I4.32p may become negative and we would 
then loose the positivity of the coefficients of the gradients squared. Nonetheless we can 
use the “strong convexity” of the one body term T{ui) in (14.261) when Ui is away from 0 
to prove: 

Theorem 10 Let h^xt £ let A he small enough, then the inf of is equal 

to the min of over homogeneous h. 

Theorem 1101 is proved in Appendix [J] 


4.5 Proof of Theorem [T] 


Using Theorem 0 and Theorem 1 101 we will next prove Theorem [TJ We thus know that 

,2 


lim lim --n-logZj’U . = lim max [- \-hext'nT'— eim)] 


where, letting A = I x 


, / \ -L , V—^ (rr^ 

TaT ^ a{x,i)=im^ for aWiGl'^ ^ 


Denote by 


/A,£M = -^log J] lE(.,)eA 


cr(fc,2)=| A|m^ 




(4.36) 


(4.37) 


(4.38) 


^€{- 1 , 1 }^ 

the finite volume free energy of the system with only vertical interactions. Thus lim^^oo f\,ii'>^) 
f\{m). We obviously have —'ipx/im) < —fx/{m) and by classical results on the Ising model 


- /A,K"i) < -f\{rn) + - 


so that 


1 mf 

lim sup lim -logZ^T r < 1™ max [-h hextm — fxifn)] 

.^^0 |A| - l^oom&Mi-.\m\<m+^ 2 ^ 

which proves that the left hand side of (12.51) is bounded by its right hand side. 
We are thus left with the proof of the reverse inequality. Let 

rh = arg minj - hextm - ^ + /a("i)| 


(4.39) 


(4.40) 


= max |m G Aii : m < m| 




h=hW 


= 


(4.41) 


(4.42) 


where Px,i{h) = ^ log Zx,h,i and Zx,h,e is the partition function of the Ising model in 

[1,£] with n.n. interaction of strength A and magnetic field /i; px{h) is the corresponding 
pressure in the thermodynamic limit i ^ oo. It is well known that 

c 


sup|pA,£(^) -Pxih)\ < . 
h ^ 


(4.43) 


13 







Then 


+ 






cr{x,i)=im(-^) ^ 




} (4.44) 


is a lower bound in the asymptotic sense for log ^ By the equivalence of ensem¬ 
bles, see Theorem [6] in Subsection IF.21 the lower bound becomes 


+ 4 - 1 - +px£ 


(4.45) 


where by (I4.43|] we can also replace by p\{h^^y. By compactness converges 

by subsequences and if £k is a convergent subsequence there is h so that —)• h and 

consequently 

=Pa(/i) (4.46) 

k 

In fact in general any limit point of p'xi{h^^^y lies in the interval ^^PaC^)] of 

its left and right derivatives, but since we are considering a d = 1 system such derivatives 
are equal to each other. Moreover by the choice of 


Hence if —)■ h then 


p'x/{h^^'^) = ^ di. 


P\{h) = m 


(4.47) 


(4.48) 


and since there is a unique h such that p'x{h) = m it follows that for any convergent 
subsequence —)■ h and therefore —?■ h. Thus the expression (j4.45p converges to 


jdexthi 


+ 


(rhy 


— hrh + px 


(4.49) 


which concludes the proof because —hrh + px{h) = fxirh). 


5 Monomials are sum of gradients 

In this section we prove a combinatorial lemma. Theorem [TT] below, which is the key 
ingredient in the proof of the gradient structure of the hamiltonian. The whole section is 
self contained and can be read independently of the rest of the paper. 

Theorem 11 Let u = (ui,.., Uk) G n = (ni, ..,nk) € and 

k 

Mn{u) = u!y ■■■vyy, y^^ni=: N 

i=l 

a monomial of degree N in the k real variables ui, ..,Uk- Then for any N >2 

k 

Mnju) = Piuf + dij{u){ui - Ujf (5.1) 
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where {pi, ■■,Pk) is a probability veetor, its component pi depending on n; dij{u) are poly¬ 
nomials of degree N—2 with negative coefficients which depend onn and there is a constant 
c so that for any positive U <1 

sup \di^j{u)\ < cU^~'^N^ (5-2) 

\ui\<U,i=l,..,k 


Proof Call Nj = ni + • • • + Uj, j = 1, k, so that = N. We will prove the theorem 
with 

N^-l 

dij{u) = (5.3) 

m=l 

with coefficients Cij-m, i<j,^<'m<Nj — 1, which depend on ni, ..,nj and satisfy the 
bound 

\ci,r,m\ < CNf (5.4) 

(j5.2l) follows from (I5.3I) - (I5.4I) which also show that the representation (|5.ip of Mn{u) is 
not unique as we can commute the factors in the monomial Mn{u) without changing 
its value. 

The proof of (|5.ip generalizes the equality 


uv 


1 


- 1 - (u — v) 

2 2 2 ^ ' 


2 


In fact we use the above identity to rewrite the factor uiU2 in M„ = Mn{u) getting 

111 2 

Mn = 2■^>21+61—62 2”^— (^1 ^2) (^•^) 

where 

hi = (1,0, ..,0), 62 = (0,1,0, ..,0) 

We thus have 

11 1 2 

2 -^^> 1 + 61—62 2 ”^— ^dn — l^dVln—e^—e 2 ("^1 ^ 2 ) 

which reminds of the discrete version of the equation Au = / that we will solve by 
iteration. There is a nice probabilistic interpretation under which the terms ^M„_|_e^_e 2 
and will be interpreted by saying that ni —)■ ni ± 1 with probability 1 / 2 , see 

the process n{t) defined below. With this in mind we introduce a Markov chain f,{t),t € N, 
f{t) £ n, where 


fc-i 

11 = Hi, Qi = {(i, x) : 1 < x < W+i ~ 1 }, i < k — 1 
i=l 

Hfc-i = {{k — l,x) : 0 < X < Nk} (5-6) 

The transition probabilities P{-,-) are all set equal to 0 except 

P{ii,x), {i,y)) = ^, \x-y\ = 1 ; 

P((i, 1), (i + 1, W+i)) = P((i, W+i - 1), (i + 1, W+i)) = i 

P{{k - 1,0), {k - 1,0)) = P{{k - 1, Nk), {k - 1, Nk)) = 1 (5.7) 
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The first line in (15.71) describes the motion on the components Oj of the second one the 
jump from to flj+i (the reverse jump having 0 probability) while the last line says that 
the endpoints of i^k-i are “traps”, namely once the chain reaches those points it is stuck 
there forever. 

We start the chain at time 0 from 


m = ihni) 


(5.8) 


We call Ti, i = 0, ..,k — 2, the first time t when ^(t) S (tq = 0) and define for i > 1, 
iTj = ± if ^(r* — 1) = A^i+i, respectively ^(xj — 1) = 1. For Ti < t < Xj+i ^{t) is a simple 
symmetric random walk, thus, by classical estimates, there are constants 6 > 0, c > 0 so 


that 


P[Ti+i -Ti> s]<ce 




(5.9) 


To establish a connection with Mn{u) and the iterates of (j5.5[) we introduce new processes 
n{t), f{t) and a{t) which are all “adapted to the canonical filtration calling 6 {t) 
adapted to the canonical filtration T) if 6 {t) is determined by {>^( 5 ), s <t}. 

Let n{t) = (ni(t), ..,nk{t)) be defined as follows. When t < n we set 


ll{i) = {C{t),N2 - ?(t), ns,.., nfc) 


For t = Ti: 


n('ri) 


In the interval ri < t < r 2 


(iV 2 , 0 ,n 3 ,..,nfc) if (Ti = + 
( 0 ,fV 2 ,n 3 ,..,nfc) if fji =- 


^ J 0 , N 3 - C(t), n4,.., nfc) if (Ti = + 

~ - ^(t),n4..,nk) if c7i =- 


(5.10) 


(5.11) 


(5.12) 


By iteration the definition is extended to all t S N. The process n{t) is indeed quite simple: 
fix 2 < z < A:, then ni{t) = ni{0) for t < Ti -2 after that it performs a simple symmetric 
random walk with absorption at 0. In the time interval Ti -2 < t < Ti-i all nj{t) = 0 with 
j < i except one, whose label is denoted £i, which jumps with opposite sign as ni{t). For 
j > i, rijit) = nj( 0 ) in Xj-i <t <Ti. 


The process f{t) is defined as 


fit) = ]Jn, 

i=l 


ni{t) 


(5.13) 


while a(t) is defined by setting 


a{t) 


1 r Ni-ne {t)-l 


j>i 


Ti-2 <t <Ti-i 


(5.14) 


We are going to prove that /(O) = M„(u) = is equal to 


t-i 


s=-l 


/(O) = E fit) ’ «(-l) = 0, t > 0 


(5.15) 
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(j5.15p will be proved by showing that 


t-i 


■= f{t) - E a{s) is a J^t-™artingale 


s=-i 


namely that 

E[m{t + 1) I Et] = 0, E[f{t + 1) - fit) I Et] = ait) 

Since we are conditioning on Et we know the process till time t, suppose that rj _2 < t < 
rj_i, call nit) = (n'^,ni+ 2 ,n^), so that all n'- = 0 with j < i except ii, while 
n'j = njiO) = Uj for j > i. Then, by (15.131) . 




u- 


j=i+l 


and 


E[fit+l)\Et] =-{u^: u.* u-^ } |J_ 




j=i+l 


SO that E[fit + 1) I Et] — fit) is equal to 

1 9 9 

^{Uii+Ui u 




Ue U- ' 


U, 


j = i+l i=*+i 

which is equal to a(t). Thus E[fit + 1) — a(t) | Et] = 0 and therefore m(t) is a martingale. 
Since P[Tk-i < oo] = 1 we can take the limit as t ^ oo in ()5.15l) which yields ()5.ip with 


1 1 


(5.16) 




= =-En =+,i < n < j - 1 

X = m, x°(s) e [1, A(j'— 1], s < t 


t>o 


(5.17) 


where x^(t) is a simple symmetric random walk which starts from Nj_ 


j-i- 


A Proof of Theorem [2] 

We preliminary observe that for any /text > 0 there is m so that hext + m = f'^im): in 
fact hext + rn — f'^ini) is positive at m = 0 and negative as m —>• 1 with /)^(m) continuous. 
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If there are several m for which the equality holds we arbitrarily fix one of them that we 
denote by we shall see a posteriori that there is uniqueness. To compute the left 

hand side of (j3.2p we introduce an interpolating hamiltonian. For t S [0,1] we set: 

+ (A.l) 

(x,i)eA 


Denote by the partition function with hamiltonian H^, by Pt,'y,L the Gibbs measure 
with hamiltonian and by Et,'y,L its expectation, then 


f 


- log 


7:^e 


(A.2) 


The thermodynamic limit of logZ£/|A| is the pressure of the d = 1 Ising model with only 
vertical interactions and magnetic field /iext + ^hext > thus, by the choice of '■ 


lim 
L —^oo 


logZ° 

|A| 


(/text + mh 

ext )mh 

ext - fx{mh ext } 


(A.3) 


To compute the left hand side of (|3.2I) we need to control the expectation on the right 
hand side of (IA.2D that we will do by exploiting the assumptions on /text which imply 
the validity of the Dobrushin uniqueness criterion as we are going to show. The criterion 
involves the Vaserstein distance of the conditional probabilities i) | {o'{y,j)}] 

of a spin a{x,i) under different values of the conditioning spins j), {y, j) / (a:,*)}. 
In the case of Ising spins such Vaserstein distance is simply equal to the absolute value of 
the difference of the conditional expectations and the criterion requires that for any pair 
of spin configurations outside (x, i) 


\Et,^,L[(T{x,i) 1 {cT{y,j)}] - Et^^^L[o'{x,i) \ {cT'{y,j)}]\ 

- '^r{x,i;y,j) <r < l (A.4) 

y,j y,j 

Since 

Et,^,L[crix,i) \{cr{y,j)]] = tanh |t ^ y)(T(y, i) + (1 - 

y^x 

+A[fj(x,i + 1) + cr{x,i - 1)] + /iextj 

{J'y,Lix, y) the kernel J^(x, y) with periodic boundary conditions in A) one can easily check 
that ()A.4p is satisfied with r as in (13.11) and r{x,i;y,j) = r^^i{x,i;y,j) with 

^7,L(a'D) i/j j) — cosh (/iext 1 2A) (x, y)lj=j + ^^x=y;j=i±lmodL^ (A.5) 


By the Dobrushin uniqueness theorem there is a unique DLR measure Pt^y which is the 
weak limit of Pt,y,L as L —oo. We denote by rrit^-y^L and rrixy the average of a spin under 
Pxy,x and P^y. We call and the measures Pt,y,L and P^y when t = 0, thus is 
the Gibbs measure for the Ising system in A with hamiltonian and magnetic field 

/iext + iii/iext) denoting its thermodynamic limit. We then have 

lim = "it,7, lim mo,7,L = (A.6) 

L—>-oo L—)-oo 
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It also follows from the Dobrushin theory that under Pt,'y,L the spins are weakly correlated; 
lei z ^ X then 


'n.>Oyi,j\,..,yn,jn 

■ ■■ry,L{yn-l,jn-l]ynJn)r^,L{ynJn;Z,i) (A.7) 

where the *sum means that all the pairs {yk^jk)^k = 1, ..,re must differ from {z,i). Thus 
there is a constant c so that 

\Et^^^L[{cr{x,i) -mt,.y,L) \ criz,i)]\ < cy (A.8) 

and also (after using Chebitchev) 


E, 


t,-y,L 


[l X] y)icr{y, i) - mt,^,L)\ 

y 


< cy 


(A.9) 


We can also use the Dobrushin technique to estimate the Vaserstein distance between Pt^-y^L 
and A- The key bound is again the Vaserstein distance between single spin conditional 
expectations. We have 

Et,^,LWi^P) I WiyJ)}] - E^o^[cr{x,i) \ {a'{y,j)}] 

= tanh ^t'^J^^L{x,y)a{y,i) + (1 - + X[a{x,i + 1) + a{x,i- 1)] + hextj 

yytx 

- tanh I + \[a'{x, i + 1) + a'{x, i - 1)] + hext + "ihext} (A.10) 

thus, calling A := cosh“^(hext — 1 — 2A), we can bound the absolute value of the left hand 
side of (lA.lOl) by: 

^ r^^L{x,i]x,j)\a{x,j) -a'{x,j)\ + At\'^J^^L{.x,y)a{yA) - rrihA 

j=i±l y 

After adding and subtracting to each cr{y,i) and recalling that Yly A,Lix,y) = 1, 

we use the Dobrushin analysis to claim that there exists a joint representation Vt,'y,L of 
Pt,-y,L and A such that 

- (x'{x,i)\] < - o-'{x,j)\] 

j=i±l 

+At(£t^^^L[\ J^ix,y){a{y,i) - m-t,7,L)l] + - mh,j) (A-ll) 

y 

Since '^y J^{x, y){cr{y, i) — mt^-y^L) does not depends on a' we can replace the £t,'y,L expec¬ 
tation by the Et^^^L expectation and after using (jA.9p we get by iteration 

Aij / \ 

St,^^L[\(y{x,i) -a'{x,i)\] < Yy^\cy+ \mt,^^L-mhA) (^- 12 ) 

with r as in ()3.ip . Since — "Io,7,l1 < ^t,'y,LAi^p)\ ~ ^'{xA)\]-, dA]) yields 

- ?t^o,7,l| < + At,l,L - mo,7,L| + |mo,7,L - rUh^^S) 
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(A.13) 


By (IMI) so that 

- mo,7,L| < ^(c7 +l"io,7,i - m-Ziextl) 

\'mt,^,L - mh,J < - mo, 7 ,L| + (cy + |mo, 7 ,L - mh,J) 

Thus rrit^^^L ^/lext ^s hrst L —)• oo and then 7 —> 0. This holds for all t and in particular 
for t = 1 hence properties (i) and (ii) are proved. Moreover, since converges 

as 7 —)• 0 to the latter is uniquely determined, as a consequence the equation 

hext + rn = fxim) has a unique solution rrih^^^ which is the limit of as 7 —>• 0. To prove 
(iii) we go back to ()A.2p and observe that 

(ai,2)GA 


Therefore 



(a:,2)GA 


(x,2)EA 

y^x 


- < |A|c 7 


(IA.2P and (|A.3[) then yield (13.2p because rtit^-y^L —t rrih^^^ as L —)• oo and then 7 —)■ 0. This 
is the same as taking the inf over all m because we have already seen that /text= /a(^) 
has a unique solution. 


B Proof of Theorem [3] 


Following Lebowitz and Penrose we do coarse graining on a scale £, i the integer part of 
7 “ 1 / 2 _ Without loss of generality we restrict L in (12.41) to be an integer multiple of i. We 
then split each horizontal line in A into Lji consecutive intervals of length I and call X 
the collection of all such intervals in A. Thus 


is the set of all possible values of the empirical spin magnetization in an interval / G X. 
We denote by M the set of all functions rn = {m{x, i), (x, i) G A} on A with values in 
which are constant on each one of the intervals I of I. Due to the smoothness assumption 
on the Kac potential there is c so that for all a, 7 and L 


(x,i)GA 


— m X, i\a 


)) 


< c7^/^A 


(B.l) 


where, denoting by Ix,i the interval in X which contains (x,i). 


m(x, i\(T) 




(B.2) 
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Thus m{x,i\a) does not change when {x,i) varies in an interval of I and therefore m = 
{m{x,i\a), {x,i) € A} S M. Then the partition function 


^7,l 


E 

m£M 


e 2 


5 Ei.a;,!, ■^i,Lix,y)m(x,i)m(y,i)+J2^ . he^tm(x,i) 


El 


m(-\a)=m 




(B.3) 


has the same asymptotics as ^ in the sense that 

1 


lim lim , . , 
7—>-0 L —>00 I A| 


log - log 


= 0 


(B.4) 


We next change the vertical interaction by replacing 

—Xa{x, ni)a{x^ + 1) —)• —Xa{x, n£)a{x, (n — l)i + 1) 


and call the new vertical energy. We then split each vertical column into intervals 

of length calling I' such intervals and A the squares I x I'. Let A = I x I', the 
restriction of m to A, so that m^{x, i), x G I, i G I' is only a function of i with values in 
Aii- Recalling the definition (j4.3[) of (t)i{m/\) we have that Z^^l has the same asymptotics 
as 







(B.5) 


meM 


where Ax^i denotes the square A which contains {x,i). 

The cardinality of M is ^ hence has the same asymptotics as 


/vmax 

^7,L,^ 


m£M 


(B. 6 ) 


Recalling the definition (j4.5p of Z^^, we are going to show that 

TlogZ«“, = Ti„gZ»» (B7) 

To prove ()B.7h we write 



m{y, 



2 


and use that J^^Lix,y) = 1. In this way the exponent in the right hand side of (IB.OP 
becomes a sum over all the squares A of terms which depend on plus an interaction 
given by 





i,x,y 


Due to the minus sign the maximizer is obtained when all tua are equal to each other and 
to the maximizer in (|4.5I) . To complete the proof of (|B.7p we still need to prove the bound 
on the magnetization: 


Proposition 2 There are Aq > 0 and m+ < 1 so that for any A < Aq the maximum in 
(jB.Gp is aehieved on eonfigurations such that for all {x,i) G A, \mAix,i)\ < m+. 
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Proof Given /i > 0 let S{m) be the entropy defined in (12.7p and let m/j be such that 


- [S'{mh) + ruh] = h 


(B.8) 


Call m* the value of rrih at h*, h* as in ()4.ip and choose m+ > m*. Fix any horizontal 
line i in A, take a magnetization m* such that m* > m+, it is then sufficient to prove that 
for all cj(x, i + 1 ) + cj(x, z — 1 ) =: hi{x), 


^-eu{mi) 


El 


a-(x)=£mi^ 




El 


E cr(x)=£m* 




(B.9) 


m 


where U{m) = — - - hextm. Since \hi\ < 2, this is implied (for £ large enough) by 


- U{mi) + S{mi) + 4A < -U{m*) + 5(m*) 
Since > m* and hext < h*, (IB.IOD is implied by 

+ h*mi + S{mi) + 4A < —^-h h*m* + S{m*) 


(B.IO) 


(B.ll) 


The function m? + S{m) + h*m is strictly concave in a neighborhood of m* where it reaches 
its maximum, hence (recalling that m* > > m* 

+ 1 ^*^* + - (^ + 
is strictly positive and (jB.9p follows for A small enough. 


C Cluster expansion 

In this appendix we will study the partition function defined in (14.8p using cluster 
expansion. 

C.l Reduction to a gas of polymers 

We shall first prove in Proposition [3] below that ^ can be written as the partition 
function of a gas of polymers F. The definition of polymers and the main notation of this 
section are given below. 

• F = {C,S,X) denotes a polymer, C its spatial support, X and S its specifications. 
C is a collection of pairs of consecutive points {L and 1 being consecutive points), 
and calling connected two pairs if they have a common point, then C is connected. 
We write i £ C or sometimes z £ F if z is in one of the pairs of C. Each pair in C is 
either a X-pair or a 5-pair, 5 and X* are the collection of all the 5 and respectively 
all the X pairs. X is the set of all points z which belong to one and only one of the 
X-pairs. 
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• \C\ is the number of pairs in C, ISI the number of pairs in S and |X| the number of 
points in X. It follows directly from its definition that |X| is even. 

• r and r' are compatible, F ~ F', if the spatial supports of F and F' do not have any 
point in common. 

• rc(F) is the weight of the polymer F. If C consists of all the possible pairs (so that 
\C\ = i) and S' = 0 then we set 

w{T) = sinh(A)^ (C.l) 

Otherwise: 

w{C, S, X) = sinh(A)l‘^l ^ = tanh{hx) (C.2) 

Each X-pair in F contributes to the weight of F by a factor sinh(A) while each S-pair 
contributes with a factor [cosh(A) — 1], as it readily follows from (jC.2p . The dependence 
of the weight on hi is through the terms Ui, i £ X. 


Proposition 3 Let F and rc(F) be as above, then 

=E n (c-3) 

£ rer 

where the sum is over all collections F = Fi,.., F„ of mutually compatible polymers. 


Proof We write 

^hcji 

^4 = Efll -1 + 11} 

(7 i i 

By expanding the last product we get a sum of terms each one being characterized by the 
pairs {i,i + 1) with — i. We fix one of these terms and perform the sum over a. 

We call cluster a maximal connected set of pairs with — 1], this will be the spatial 

support of a polymer. The sum over u factorizes over the clusters. After writing 

gAo-io-i+i _ ^ _ sinh(A)(TiC7i+i + [cosh(A) — 1] 

we call (i,i + 1) a X-pair if it has the term sinh(A)(Ti(7i+i and a S'-pair if it has the term 
[cosh(A) — 1]. Notice that if i belongs to two X-pairs then we have a product of two Uj 
which is equal to 1. Thus the sum over the spins in a cluster C becomes a sum over tc(F) 
with the spatial support of F equal to C. In this way we get (IC.3p . 


We shall also consider the partition function 

^; = En^i(r) (C.4) 

£ rer 

where tci(F) is obtained from rc(F) by putting Ui = 1. 
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C.2 The K-P condition 


The K-P condition for cluster expansion requires that after introducing a weight |r| then 
for any T 

|u;(r')|el^'l < |r| 

rvr 


Proposition 4 For A small enough we have that 


^ < |r|, 6>0 

r'/r 


with 


|r| = |C'(r)|+ 1 , e^=A-5/^2 

having called C'(r) the spatial support ofT. 


(C.5) 


(C.6) 


Proof We are first going to prove that for A small enough 

Y < 1 

r':C'9l 


(C.7) 


Fix C and consider all F with spatial support C, i.e. C'(F) = C, so that |F| = IC*! -|- 1 =: n. 
Then 


^ u;i(F)e(^+'’)l^l < e(^+^)[sinh(A)e(^+^)]"-^(l + 

r:C(r)=c 


-if, , [c osh(A) - 
sinh(A) 


(C.8) 


Therefore the left hand side of (IC.7P is bounded by 


n>2 


which vanishes when A —)• 0, because by (IC.6P Ae^^ vanishes as A —0. Hence (|C.7h holds 
for A small enough. 

To prove (|C.5p we first write 


E I^(r0|e 


{i+b)\r'\ 


< 


Y 




(C.9) 


rvr 


rvr 

and then use (IC.7p to get 

rvr iGC(r)r':C(r')3i 


E < E E < |r| 
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C.3 The basic theorem of cluster expansion 

The theory of cluster expansion states that if the K-P condition is satisfied then the log of 
the partition function can be written as an absolutely convergent series over “clusters” of 
polymers. To define the clusters it is convenient to regard the space {T} of all polymers as a 
graph where two polymers are connected if they are incompatible, as defined in Subsection 
1C.II Then a clnster is a connected set in {T} whose elements may also have mnltiplicity 
larger than 1. We thus introduce functions I : {T} —)■ N such that {T : /(T) > 0} is a 
non empty connected set which is the cluster defined above, /(T) being the mnltiplicity of 
appearance of T in the cluster. With such notation the theory says that 


logZl^ = YW^^ 

I r 

(C.IO) 

log Zf = Ywl, Wl-.= 

(C.ll) 


I r 


where the snms in (IC.lOp - dC.llI) are absolntely convergent. The coefficients o/ are com¬ 
binatorial (signed) factors, in particular o/ = 1 if I is supported by a single T. We will 
not need the explicit expression of the aj and only use the bound provided by Theorem 
[12] below. We nse the notation: 

|/|i = j;/(r), ||/|| = ^|r|/(r) (c.12) 

r r 


Theorem 12 (Cluster expansion) Let A be so small that the K-P condition (1C.51) 
holds. Let T be a polymer and I a subset in {/} such that /(T) > 1 for all I G X (I 
could be the whole {I}). Then 

supe-^H-^H (C.13) 


Observe that the absolute convergence of the sum in (IC.lOD - dC.lip is implied by ()C.13I) 
with I = {I : I{T) > 1} as it becomes 

Y (C.14) 

7:/(r)>l 

becanse inf/gj as the inf is realized by I* which has /*(r) = 1 and = 0 

for all T' / r. (jC.lip proves that the sum in dC.llI) and hence the sum in (|C.10p are both 
absolutely convergent. 


D Proof of Theorem |4] 

In this section we will prove Theorem [5| as a direct consequence of Theorem [T21 
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D.l Proof of (l4.10p 

We start from (jC.lOp and observe that 


:= ai = {ai JJ 

r r r 


The last factor is 


equal to where N{-) is determined by I: 


N{x) = Ewi x&x{r) (D.15) 

r 

hence p4.10l) . Recalling (I4.12p we observe that |A^(-)| is even because the cardinality of 
each X(r) is even. 


D.2 The term with |A^( )| = 0 

The term with |A^(-)| = 0 is a constant Aq (i.e. it does not depends on u) and it will not 
play any meaningful role. It is bounded as follows: 


Lemma 2 There is a constant c (independent of u and i) such that 

l^ol < 

Proof By ficm 


i-^"! s EE E I*"'! 

i=l CBi I J(C,C,&)>0 

< y^y^[cosh(A) — < fcA^ 

i=l CBi 


D.3 Proof of fl4.15p 

We have 

o-j-iUiUj = lx={ij} W 

T={C,S,X) 7:/(r)=l;7(r')=0 if X' ^ 0 and V 

Thus by (IC.13P 

\aj-i\ < ^ lx={ij}Wi{T)e\^\ 

T={C,S,X) 

< lx={ij}(sinh(A))l*“^l(cosh(A) - 

r=(c,s,x) 
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which is bounded by 


\aj-i\< Y. (smh(A))l*-^l(cosh(A) - 

n>0,m>0 

We have thus proved the second inequality in (I4.15p . 

To prove the first one we call T* = (C*, (S'*, X*) where C* = {i, i+1), S* = 0, X* = {i, i+1} 
and write 

Oi^i+iUiUi+i = sinh(A)ttirii+i + '^ 

/:|/|l>l;/(r*)=l;7(r')=0 if T' ^ T* and 

+ E U.(.,n E (D-“) 

r=(c,s,x),r^r* /:/(r)=i;7(r')=o if x' 0 and r' ^ r 

If I is the second term on the right hand side then 11/| | > 2 + 2 so that this term is bounded 
by 

|rc(r*)|el'" < sinh(A)e^e“^^ 

Proceeding as in the proof of the second inequality in (14.151) we can bound the last term 
on the right hand side of (ID.lbj) by 

< ^ (sinh(A))(cOsh(A) - l)n+m^ 2 +n+m+l 

n>0,m>0,n+m>0 


which proves the first inequality in ()4.15p . 


D.4 Proof of (I4.13P 


If I determines N{-) then for all j 



NU)< 

r:C(r)9i 

(D.17) 

hence 

R{N{-))< |r| 

r:/(r)>o 

i i r:C(r)9i 

(D.18) 

Thus 


l|/||>l|A^(-)ll 

(D.19) 

so that the left hand side of (I4.13P 

is bounded by: 



E E 

|1P/| < 

(D.20) 


rgi/:/(r)>o,||/||>M rsi 


having used (IC.lSp . (I4.13P then follows from (1C.71) . 


27 




















E A priori bounds 

We will extensively use the bounds in this section which are corollaries of Theorem HI 


Corollary 1 There are eonstants Ck, k > 0, so that for any i S {1,..,^}, k > 0 and 
M > 4, 

^ IA(-)f I^Ar(.)l < CkM^e-^^ = CfcA5/3e-^(^-4) (E.l) 

N(-):N(i)>0,\N(-)\>2,\\N {■)]]> M 

Proof 

It follows from Theorem HI see (j4.13p . 


Corollary 2 There are constants c'j^, k > 1, so that for any i and i G [1,^] 


E 


d‘ 


\k—l 


d 


dui^ ■ ■ ■ duij^_^ duj 


logZ;j<4A 


for any A as small as required in Theorem\^ Moreover 

Ti(tt) = 0 if \ui\ = 1 


(E.2) 


(E.3) 


Proof We write logZff^ = Ki + K 2 where Ki is obtained by restricting the sum on 
the right hand side of (|ClOp to |1V(-)| ^ 2, K 2 is the sum of the remaining terms. By 
(j4.14p - (l4.15p we easily check that Ki satishes the bound in (IE.2p . We bound 


E 


dui^ ■ ■ ■ 


A 

dui 


K 2 


by 


A E \N{-tR{N{-)Y 

M>2 ||Ar(.)||=M,Ar(i)>0 


(jE.2p then follows from (jE.ip . ()E.3p follows directly from the definition of Tj(rt). 


Corollary 3 Recalling (j4.14p and writing a = oij-i, 
A ctj-iUiUj = « E '“A ^ E 

i<j i j>i 


(E.4) 
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F Proof of Theorem [5] and Theorem [6] 

We write ||?;|| for the sup norm of the vector v: ||u|| := \vi 


F.l Proof of Theorem [5] 


Existence. By (IE.2h we can use the implicit function theorem to claim existence of a small 
enough time T > 0 such that the equation 


m = u{t) + t'^{u{t)) 


(F.l) 


has a solution u{t),t E [0,T], such that: u(0) = m, u(t) is differentiable and ||tt(t)|| < 1, 
recall that ||m|| < 1. 

If A is small enough ()E.2p with k = 1 yields 


d 


max sup It:—'I' j(u)| =: r < 1 


(F.2) 


so that the matrix 1 + tV'I'(u(t)), (V'I')ij = is invertible for t < min{T, 1} and 

therefore for t < min{T, 1} 

t((t) =/(n(t), t) := — ^1 + tV'I'(u(t))^ 'h(u(t)), u(0) = m (F-3) 


By (IF.2p - ()E.2l) f{u,t) is bounded and differentiable for t < 1 and ||ri|| < 1, thus we can 
extend u{t) till min{l,T} where r is the largest time < 1 such that ||u(t)|| < 1 for t <t. 
Thus for t < T (jF.ip has a solution u{t) which we claim to satisfy ||u(t)|| < 1. To prove the 
claim we suppose by contradiction that there is a time t < r and i so that \ui{t)\ = 1. By 
dET]), rui = Ui + t'^i{u) = Ui (having used ()E.3M . We have thus reached a contradiction 
because ||m|| < 1. Thus the claim is proved and as a consequence r = 1 and therefore we 
have a solution of (jF.ip for all t < 1 with. 

Uniqueness. Suppose there are two solutions u and v. Then 


u — u = T(u) — T(rt) 

Define u{s) = su + (1 — s)u, s € [0,1], then 

ll^t —'y||< [ ||VT(rt(s))(u — u)|| ds 

Jo 

Since ||u(s)|| < 1 by ()F.2p ||VT(tt(s))(n — u)!! < r||ti — u||, so that ||n — u|| < r||ti — u|| and 
therefore u = v. 

Boundedness. Calling u = u{t) when t = 1, by (jF.ll) and (|E.2p 

11^*11 < ||"i|| + ||'l'(u)|| < ||m|| +ciA (F.4) 

so that if ||m|| < m+ then for A small enough ||rt|| < 1 and therefore there exists such 
that ||h|| < /i+. 
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F.2 Proof of Theorem [5] 

Since 


1 


- (^£(m) = ^log{e lm(»=m 

<7e{-i,i}^ 


p- J2x i{->^o-{x,i)aA(x,i+l)-hia{x,i)} 

=m^ ’ 


we have for free 


1 A,, (m) 




} (F.5) 


(F.6) 


and we are thus left with the proof of a lower bound for —(j)p(m). 
Call li = {(x, i) : X < £ — let a' € (^,a) and 

Bi = {a{-,i):\ Y - mi]\ < £'"'} 

(x,i)Gli 


(F.7) 


Let /i be the Gibbs probability for the system with vertical interactions and magnetic 
fields h. We look for a lower bound for 


/i [^{Pl Bi} n {m(-|cr) = m} 


By the central limit theorem 




B^ 


< e 


-bi^ 


, b>0 


because the spins in Jj are i.i.d. with mean m*. Moreover 




{m{-\a) = m} | 


> e 


-4A£l + “r,-tl + “ 


(F.8) 


(F.9) 


because, given {piSj}, there is at least one configuration in the complement of li on each 


horizontal line. Thus 


n M-k) = m]\ > (1 - 


hence 


-Urn) > - llog{(l 


which together with (IF.6P proves (I4.18p . 


G Proof of Lemma [T] 


We first write 




efr 




f- 2 


Yi-^ - (^ext - h,)u^ - \og{e^^ + 


2 = 1 




+ ~ Ui — hextl^i — log + Aijj 

2=1 


(G.l) 
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We have log(e^* + e = hiUi + S{ui), the entropy S{u) being defined in (I2.7I) - (I2.8I) . 
Thus 

^ \I /2 

^ 1,1 = 

2=1 

The term with /lext^i in ()G.2p becomes 

^ext ^ ^ ^2 T -^^ext ^ ^ ^ 2+1 T 2A/iext ^ ^ '^i 

i i i 

which can be written as 

- hext E + Ui+i)(Ui+i - (G.3) 

2 2 

After an analogous procedure for the term with {hi — Ui)^i we get (|4.26p . 


H Proof of Theorem [8] 

We say that a function T(u) is “sum of one body and gradients squared terms” if 

e 

F'iu) = E + E - '^jf 

i=l l<i<j<^ 

for some functions f{u) and bij{u). Thus (I4.28|) claims that is “sum of one body and 
gradients squared terms”. We say in short that the “gradients squared terms are bounded 
as desired” if 

E ^ cA^+i E('“* “ "“m)^ 

l<i<j<£ i 

Hence (j4.29p will follow by showing that the gradients squared terms of are bounded 
as desired. 

We will examine separately the various terms which contribute to and prove that 
each one of them is sum of one body and gradients squared terms and that the latter are 
bounded as desired. 


H.l The 0 term 

By (031) 

0 = E ^ '^{Ui+l - Uif 

Af (')^0 i=l 

Call 0^^) the above expression when we restrict the sum to N{-) : |A^(-)| = 2 and call 
0 (> 2 ) = 0 _ 0 ( 2 )_ Thus is equal to the sum of Ajv(.) over N{-) : |A^(-)| > 2, i.e. 
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|-^(•)l — recall in fact from Theorem 0] that ^Ar(.) = 0 if N(-) is odd. We start from 
0(2) which, recalling (IE.4p . is equal to 

0(2) = Q, ^ ^2 _ 1 - X){ui+i - Uif - i - mf (H.l) 

i i 

Thus 0(2) is sum of one body and gradients squared terms. To prove that the latter are 
bounded as desired we write 

{uj - Uif < {j - i) ^(ufc+i - Ukf (H.2) 

k=i 

and call n = k — i > 0, m = j — k > 1. We then use (14.151) to bound the sum of the terms 
with the gradients by 

J](ufc+i-Ufc)2{cAe-2'>+ ^ (m + n)cA”^+’^e"^+"} (H.3) 

k n>0,m>l,m+n>l 

which is the desired bound because | < §• 

We rewrite 0 (> 2 ) using (j5.ip for each one of the factors Thus given N(-) we call 

ii < i 2 < ■ ■ ■ < ik the sites where N{-) > 0 and call n = (A^(ii),.., N{ik)). We then apply 
(j5.1l) with ui = Ui^,...,Uk = Ui^ so that pi and dij in (15.11) become functions of u and 
N{-). We then get 

0 (> 2 ) ^ Y p.nj^(')l + Y dij{ui-Ujf} (H.4) 

A^(-):|^(-)l>4 i:N{i)>0 j>i:N{j)>Q,N(i)>Q 

which is sum of one body and gradients squared terms. To get the desired bound on the 
latter we use (IH.2I) and (15.2p to get 

^(“<=-“'=+1)" E fO-') E c|JV(-)P|A„h|} 

k i,j-j>k>i N{-):\N{-)\>4:,N{i)>0,N{j)>0 

Since both N{i) > 0, N{j) > 0 then j — i < R{N{-)) and given R{N{-)) > k — i there are 
at most R{N{-)) possible values of j. Therefore the above expression is bounded by 

E(“*-“m)"E E ll"(■)f|.4„(,,|) 

k i<k N{-):\N{-)\>4,N{i)>0,R{>^{-))>k-i 

We upper bound the above if we extend the sum over N{-) such that 

\N{-)\ > 4:,N{i) > 0 , ||1V(-)|| > 7k-i, Ik-i '■= max{4, A: - i] 

We then apply (jE.ip with A: = 5 to get 

- Uk+if ^ ^{uk - Uk+i)^{Y^ 

k i<k k i<k 

The curly bracket is bounded by 

4^5 + '^{n + 4fe-^ < c 

n>l 

Thus also 0 (> 2 ) is bounded as desired. 
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H.2 The term /lext 
By Km 

= {l-u^)(^{ai-X){ui+i+Ui-i)+ aj-iUj+ ^ (H.5) 

j>i+l N{-)-.\N(-)\>4: 

where ei{j) = 0 if j 7 ^ i and = 1 if j = z. 

Call Qi := (1 — uf){ai — A) then the first term contributes to by 

^ - {gi - gi+i){ui - Ui+i)^ = ^ 2giUi + (ai - A) + Ui+i){ui - Ui+if 

i i i 

which is sum of one body and gradients squared terms. By (I4.15P the coefficients of the 
gradients squared are bounded in absolute value by 2 cAe“^^ which is the desired bound 
because | < |- 

By an analogous argument and writing g[ := (1 — tt?), the contribution of the second term 
in (IH.5h is 

'^aj-i{2g[ui - {g'i - g'j){ui - uj)^ = '^aj-i(2giUi + {ui + Uj){ui - 
i<j i<j 

which is sum of one body and gradients squared terms. We bound the latter using (jH.2p 
and the second inequality in ()4.15p to get 

^ 2 c(eA)'=} 

which is the desired bound because the curly bracket is bounded by c'X^. 

To write the contribution to of the last term in (IH.5p we introduce the following 

notation. Given N{-) : N{i) > 0 we call N'{-) = N{-) — e* and N''{-) = N{-) + e*. Let 
then h < i 2 < ■ ■ ■ < ik the sites j where N'{j) > 0, n = N'(ik)) and denote 

by pj, djj, the corresponding coefficients in (15. ip . Similarly let i[ < i 2 < ■ ■ ■ < i'j^ the 
sites j where N''{j) > 0, n = (W'(ii),.., N"{ik)) and denote by p'j', the corresponding 
coefficients in dSH). Then the contribution to ‘hj of the last term in (jH.Sp can be 
written as 

j: .4„,, j: iv(i)( ^ ^ 

N N{-)\> 4: r.N{i)>0 j-.N'{j)>0 j-N"(j)>0 

+ dJj,{uj-UjY- d+.,{uj-Uff}iR.6) 

j<j':N'{j)>0,N'{j')>0 j<j'.N"{j)>0,N"{f)>0 

which is sum of one body and gradients squared terms. To bound the latter we examine 
the terms with d~, those with d'^ are analogous and their analysis is omitted. For the d~ 
terms we get the bound: 

^ N{i) Y. <!|A'(-)t(%-V)" 

N{-)-\N{-)\>4 v.N{i)>0 j<j'-N'{j)>0,N'{j')>Q 

^ Y1 c\N{-)\'^{uj -Uff 

Af(.):|7V(.)|>4 j<j'.Nij)>0,Nij')>0 

which has an analogous structure as the gradient term in (IH.4P . Its analysis is similar and 
thus omitted. We have thus proved that hext Yi bas the desired structure. 
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H.3 The term V 

Z—/Z I 

We introduce the following notation: given i, N{-), N', a G {—1,1}, a' G {—1,1}, 
N{i) > 0, N'{i) > 0, we call 

iV(.) = iV(.) + N'i-), K ^ ■■= N{-) + {a + a')e. 

Then ^ • T? is equal to 

^ iV(i)iV'(i)A^(.)Gl^,(.)(-l)^+i( PjiK)uf^ 

i N{-),N'i-),a,a' r-K{j)>0 

+ Y. d3,r(K){u,,-u,f) (H.7) 

j<j':K{j)>0,K{j')>0 

which is sum of one body and gradient squared terms. Let 

Cj,3'--=Y E iV«iV'Wl^w(.)||^iV'(.)| E Kj'iK)\ 

i N{-),N'{-),a,a' j<j'-.Kij)>0,Kij')>0 

then the gradient squared terms are bounded by Ylj<j' ~ We have 

E W(i)iV'(i)|A^(.)||7l^,(.)| ^ c(|iV(.)| + |iV'(-)|+2)3 

i N{-),N'{-) j<j'.ff{j)>0,N{j')>0 

because 4 is the cardinality of (cr, ex'). Moreover 

Cj,f < ^cY Y |T^(.)||T^.(.)|(2|iV(.)|)"(2|iV'(.)|)^ 

* N{-),N'{-)-.N{j)>0,N{j')>0,N{i)>0,N>(i)>0 
By the symmetry between N{-) and A^^(-) we get with an extra factor 2: 

C,,,. < 8c 44 |T^(.)||T^q.)||iV(.)|"iV'(.)r 

i N{-),N'{-)-.N{j)>0,N{j')>0,N{{}>0,N'{i)>0 

Moreover either R{N{-)) > [j' — j)/2, or R{N'{-)) > (j' — j)/2 or both, hence 

C,,,v<8c 44( ^ l^w(.)l|W(-)l' E E 

N{-).N{j)>0,R{N{-))>i^ i:Af(0>0 N'(-y.N'(i)>0 

+ E l■^^'(-)ll'v{■)l‘ E E 

N{-y.N{j}>0 i:N{i)>0 jv'(.):Ar'(i)>0,/?(iV'(-))>T^ 

By (mi) 

< 8c 44( ^ |T^(.)||iV(.)niV(.)|c4e-2^ 

Ni-yNij)>0,R{N{-))>^^ 

+ Y l^w(.)l|W(-)nW(-)|c4e-^”-{2,Z^}^ 

N{-yN{j)>0 
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Using again (lE.ll) 


Cjj> < 8c4^2c4e-2^C5e-^““^2, Vl =; 

Hence 

E CjA^r - Ujf < Y.A+1 - uuf E (J' - 

j<j' k 3,j'-j<k<j' 

The last sum is bounded proportionally to e“'^^ (details are omitted) which gives the 
desired bound. 


H.4 The term Yhi 

Recalling (I4.27P and (I4.25P the contribution to due to is 

£ 

'^{hi-ui){l-uj){'^ aj-iUj+ E (H.8) 

i=l j>i+l N{-):N(i)>0,\N{-)\>A 

We have 

3 oo ^ CO 

(/i - m)(1 - = y - 2 ^ =■ E (H.9) 

k=2 k=l 

with \Kk\ < 1; since |u| < tt+ < 1 the series converges exponentially. We start from the 
terms with aj-i'. 

E E “j-* E = E E "j-* E '^kiipA’^^'^+PjufA + dA - u^} 

i=l j>i+l k>l i=l j>i+l k>l 

where {pi,Pj) is the probability vector introduced in Theorem 1111 and d the corresponding 
coefficient. They depend on the pair {2k + 1,1) and \d\ < ck^u^. This is sum of one body 
and squared gradients terms and we are left with bounding the latter. We have the bound 

i i 

E E A-i\'^ck*^'kkf{ui-Ujf<'^ E \aj-i\c'{ui-uA 

i=l j>i+l k>l i=l j>i+l 

which satishes the desired bound as proved in Subsection lH.il 

We next study the last term on the right hand side of (IH.8p . Proceeding as before we 
check that it is sum of one body and gradients squared terms and next prove that the 
gradients are bounded as desired. We first bound them by 

EE E N{i)\^N{-)\'^c{2k + \N{-)\fuf{uf-Ujf 

i j<j'N{-y.N(i)>0,N{j)>0,N{j')>0,\N{-)\>A k>l 

We have {2k + |A^(-)|)^ < so that we get the bound 

EE E ^(*)l^7V(.)|c'|iV(-)P(^^T 

* j<j' N{-y.N(i)>0,N{j)>0,N{j')>0,\N(-)\>A 
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with 


c' := y^^{2k)^uf 
k>l 

We can perform the sum over i to get 

j<j'N{-):N{j)>0,N{j')>0,\N{-)\>i 

We are thus reduced to the case considered in Subsection IH.li we omit the details. 


I Proof of Proposition [T] 


Recalling that £^{u) := {h{u) — n)(l — tt^), we have, supposing u' > u', 


m - ^{u) = 


du 


du < a{ui — Uj), 


(I.l) 


with a = max —. 

|u|<i du 


Thus 0i{u) < a and by (jH.9|) 


a 


max 

| m |<1 



2uY, 

k=l 


^2k+l \ 

2k + IJ 


< max 
|u|<l 



3 

8 


having retained only the term with k = 1. 


J Proof of Theorem [TO 

We shall use in the proof that in all terms but (T(u) — hextu), cf. (IG.2h . are propor¬ 
tional to A. 

Calling u the minimizer of (T{u) — hextu) : 

• It will follow from Lemma 0] that the minimizer u* of has components u* such 
that |n* — ?i| < (for all A small enough), and that the minimizer v of f{u), f{u) 
the one body term defined in (|4.28l) . is such that \v — u\ < A^/^; 

• Since the gradient of vanishes at u = (uj = u, i = 1, ■■■,(-)■, cf. (I4.28p . n is a 

critical point of 

• T{u) is a convex function and its second derivative T''{u) is a strictly increasing, 

positive function of tt G (0,1) which diverges as n ^ 1, as it follows from (j4.22p . 
Then the matrix is positive definite in the ball u ■. \ui — u\ < A^/^, cf. 

Proposition 0 
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As a consequence, the minimizer of in the ball coincides with v and since u* is in the 
ball it coincides with v, thus proving that all the components of u* are equal to each other. 
We are thus left with the proof of Lemma U] and Proposition [SJ We need a preliminary 
lemma. 


Lemma 3 For any hext S [ho,h*] there is a unique u such that 

d 


du 


{T{u) - hextu} 


= 0 


and there is ouq > 0 so that 


inf -r^T{u) 

/iext6[/io,/i*] du^ 


> Cho 


(J.l) 


(J.2) 


Proof The proof follows from the fact that the second derivative of T{u) is positive away 
from 0 and in ( 0 , 1 ) increases to oo as tt — 1 . 


Fix all Uj,j 7 ^ i and call F{ui) the energy as a function of u*. Then 

Lemma 4 There is > 0 so that for all X small enough the following holds. Let 
hext £ [^0)^*] o-iT-d u as in Lemma\^then 

inf F(ui)>F(u) + c',„Ai/2 (J.3) 

Proof By (l,T.2j) 

inf |{T('w) - hextu} - {T{u) - hextu}\ > 

We are going to show that the variation of all the other terms in (jG.2h are bounded 
proportionally to A and this will then complete the proof of the lemma. We have 

\{hi - Ui){l - uj)\ < c, (1 - u-< cA 

(the first inequality by (IH.9p . the last inequality by (jE.2p L 

Call G{ui) the value of log when tanh(/ij) = Ui and the other hj are fixed, then 

\G{ui) - G{ui)\ = \ ^ - u')| < cA|ui - u'l 

N{-y.N(i)>Q 

where, to derive the last inequality, we have used Theorem U) 


As a corollary of the above lemmas 
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Lemma 5 For A small enough the inf of is achieved in the ball u : max{|tij — iil < 


Proposition 5 For A small enough the matrix strictly positive in the ball 


u : max{|ui - \ = 1, 


Proof From Lemma [3] and Corollary [2] one obtains 




> Cho - Aci, for i = 1, 2,L 


For any i, 


from ()4.13p and Corollary [2j 






, duiduj 
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